Direct observation of correlation time of dynamic nuclear polarization in single quantum dots 
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The spin interaction between an electron and nuclei was investigated optically in a single self-assembled 
InAlAs quantum dot (QD). In spin dynamics, the correlation time of the coupled electron-nuclear spin system 
and the electron spin relaxation time play a crucial role. We examined on a positively charged exciton in a QD 
to evaluate these key time constants directly via the temporal evolution measurements of the Overhauser shift 
and the degree of circular polarization. In addition, the validity of our used spin dynamics model was discussed 
in the context of the experimentally obtained key parameters. 

PACS numbers: 73.21.La, 78.67.Hc, 71.35.Pq, 71.70.Jp 



I. INTRODUCTION 

The hyperfine interaction in semiconductor quantum dots 
(QDs) is enhanced owing to the strong 3D confinement of the 
electron wave function; consequently, this has attracted con- 
siderable attention from the fundamental and practical points 
of view. The sophisticated control of nuclear spin polarization 
(NSP) is required for fascinating applications such as a long- 
lived memory at an atomic level' and qubit conversion be- 
tween an electron spin and a photon^.. In semiconductor QDs, 
the enhanced hyperfine interaction provides the possibility of 
polarizing nuclear spins (n-spins) in one direction with the op- 
tically selective excitation of the electron spin (e-spin). In fact, 
a large NSP of up to 30-60% was observed recently in inter- 
face GaAs QDsi, self-assembled InAIAs QDs^A In(Ga)As 
QDsi— , and In? QDftiS. In these QDs, the confined electron 
is subject to a large nuclear field (Overhauser field: Bn) up to 
several Tesla. In addition, a fundamental interest is the knowl- 
edge of the decay of the e-spin polarization induced by the 
random fluctuation of the Overhauser field. This fluctuation 
induces an additional e-spin precession around the effective 
magnetic field, and it imposes an inevitable contribution to 
the e-spin relaxation and decoherence through the transverse 
and longitudinal components of fluctuation, as predicted 
in a previous worbil. 

From these points of view, it is necessary to examine the 
spin dynamics of a coupled electron-nuclei (e-n) system that 
is well isolated in a QD. In the framework of the current semi- 
classical dynamics model of NSP formation^'—, the following 
two physical quantities play a significant role; the correlation 
time of the coupled e-n spin system and the e-spin relaxation 
time. 

The correlation time Tc indicates the characteristic time dur- 
ing which the randomly fluctuating effective field is consid- 
ered to be constant in magnitude and direction according to 
the traditional spin relaxation theor y'^-'"^ , and it induces a ho- 
mogeneous broadening (2^/tc) of the target e-n levels. Since 
the NSP formation rate is very sensitive to the degree of en- 
ergy mismatch in the e-n spin flip-flop process, which is deter- 
mined by the splitting and the broadening of the correspond- 
ing e-n energy levels, has a crucial influence on the NSP 



dynamics. Despite its importance, this key quantity has been 
generally used as a fitting parameter in the model calculations 
to reproduce the observed results, and the direct experimen- 
tal estimation of has not been reported thus far in any QD 
materials. Meanwhile, the e-spin relaxation time also af- 
fects NSP formation by changing the longitudinal component 
of the e-spin polarization. As mentioned in a later section, the 
e-spin relaxation rate (1/ts) can generally be expressed as a 
Lorentzian function of the effective magnetic field, and it is 
influenced by the aforementioned correlation time through its 
width of 2fi/Tc. Accordingly, we measure the e-spin relax- 
ation time at a zero effective field Tjo, which gives the maxi- 
mum amplitude of 1 /ts; this amplitude can be evaluated inde- 
pendently of Tc. 

A secondary interest is the possibility of the measurement 
of NSP in a QD structure through the degree of circular po- 
larization (DCP) of positively charged exciton emissions. The 
DCP of the time-integrated photoluminescence (PL) has been 
used as a powerful tool to detect NSP or the Overhauser shift 
(OHS), which is the energy shift in the electronic level in- 
duced by Bn, in bulk and quantum well structures for a long 
time. However, the method to probe NSP in single QDs has 
been limited only to the change in the energy splitting of the 
PL lines; this is the simplest way to evaluated OHS, but its 
accuracy has been limited by the spectral resolution of the ex- 
perimental setup. In the coupled e-n system, by using the DCP 
of positively charged exciton emission, which directly reflects 
the e-spin polarization, it may be possible to follow not only 
the e-spin but n-spin dynamics, and the study of the coupled 
e-n system may act as a tool for sensitive measurements of 
QD-NSP 

In this study, we investigated the e-n spin dynamics in QD 
structures by using the DCP of the positively charged exci- 
ton (X^). The DCP of PL, which is mainly determined 
by the e-spin polarization, changed in synchronization with 
the OHS or the energy splitting of the e-spin levels, and this 
phenomenon provides the possibility of the sensitive probing 
of the QD-NSP. By taking advantage of this feature, the key 
quantities (Tq and Tjo) were evaluated directly from the ex- 
perimental data. In addition, we extended the dynamic model 
of NSP by including the dynamics of the X+ state, and we 
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confirmed the validity of our model by comparing the time- 
resolved OHS and DCP measurements with the calculated re- 
sults. 



II. SAMPLE AND SETUP 

Self-assembled Ino.75Alo.25As QDs, which were embedded 
with an Alo.3Gao.7As layer grown on an undoped (100) GaAs 
substrate by molecular beam epitaxy, were used in the exper- 
iments. Assuming a lens-shaped QD with a typical diameter 
of ~20 nm and a height of ~4 nm derived from atomic force 
microscopy (AFM) measurements, the number of nuclei in 
a single QD was estimated roughly to be ~3xlO'*. Micro- 
PL measurements were performed at 6 K under longitudinal 
magnetic fields (B^) of up to 5 T. A cw-Ti:sapphire laser of 
~728 nm, which provides the transition energy at the wet- 
ting layer of InAlAs QDs, was employed to illuminate the QD 
sample. The QD emission spectra were detected using a triple 
grating spectrometer and a liquid N2-cooled Si-CCD detector 
Though the energy resolution of our setup was ~ 12 jueV, it can 
be improved to 5 //eV by spectral fitting. The polarization of 
the excitation light was controlled using an electro-optic mod- 
ulator (EOM), and the Zeeman splitting energy and the DCP 
of a target single QD spectra were evaluated. 



III. RESULTS AND DISCUSSION 
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FIG. 1 : (a) PL spectra from the target single QD at zero magnetic 
field, (b) PL spectra of X+ state at 5^ =5 T with linearly (gray 
squares) and circularly (cr+/cr_ : open/solid circles) polarized excita- 
tion. The cr_(+) PL component is positioned at the higher (lower) en- 
ergy side. The difference between the Zeeman splitting AEz for the 
circularly and linearly polarized excitations is defined as the Over- 
hauser shift (AEqus), and it is evaluated as A£ohs =98 /jeV (-21 
fieW) for cr_ (cr+) excitation, (c) dependences of AEz of the 
PL line (upper panel), the energy splitting of the electronic level AE^ 
(middle panel), and the DCPs of and X" (denoted by the black 
and gray symbols in lower panel). The solid line in the upper panel 
is the calculated AEz at Sn =0. 



A. Electron and nuclear spin polarization in a coupled system 

First, we investigate the availability of DCP as a pow- 
erful measure of the electron and nuclear spin polarization in 
a QD. Fig.lifa) shows the PL spectra obtained from the target 
single QD under a zero magnetic field. Regardless of the un- 
doped sample we used, the PL spectra from the various charge 
states that originated from the same single QD were observed. 
The charge states of three peaks with high intensities were as- 
signed to a neutral biexciton (XX"), a neutral exciton (X"), 
and a positively charged exciton (X^) from the lower energy 
side, respectively. Since X^ had the strongest PL intensity in 
our QD sample and because there were no dark exciton states 
present, it is expected to be the dominant contributor in NSP 
formation. 

Fig. mb) depicts the X+ PL spectra at -5 T for the 
linearly (denoted by grey squares) and the circularly {cr+ 
and cr_: denoted by open and solid circles) polarized exci- 
tations. In the X^ states composed of the spin-paired two 
holes and an electron, the exchange interactions between the 
electron and hole spins play no role, and the energy split- 
ting of the PL lines (AEz) is determined solely by the Zee- 
man interaction of the spins with the (effective) magnetic 
fields. Under this condition, AEz can be expressed as follows: 
AEz = ^zA'bSz + gz/^B (Bz ± ^n), where g'f^^ denotes the elec- 
tron (hole) g-factor in the growth direction, fiB denotes the 
Bohr magneton, and Bn denotes the Overhauser field. Since 
the hole spin has a low probability of existence at the nucleus 



site, the effect of Bn on the hole spin could be neglected ex- 
cept for the special casai^. Under a large of a few Tesla, 
Bn manifests itself as an OHS defined as AEohs - ^z/^b^n- 
Since Bn is essentially zero for the linearly polarized excita- 
tion, OHS is deduced from the difference between AEz for 
the circularly and linearly polarized excitations, and it is eval- 
uated in Fig. [nb) as A^ohs =98 ^eV (-21 jueV) with o-_ 
(cr+) excitation. As per our definition, the cr_ (cr+) excitation 
generates Bn in the opposite (same) direction to Bz, and it in- 
duces an apparent increase (decrease) in AEz because of the 
relation gl ■ <0. Hereafter, we focus on the (t_ case where 
the compensation of Bz via Bn is achieved; consequently, the 
bistabilities of NSP have been observed for external param- 
eters such as the excitation power—, excitation polarization^, 
and external magnetic field^. 

Figure [TJc) summarizes the effects of NSP on the X^ PL 
observed in the Bz dependence measurement. In the experi- 
ment, the excitation polarization was fixed at cr_, and the ex- 
ternal field was swept from 4.0 T to 5.0 T with a sweeping 
rate of 0. 1 1 T/min. The symbols and the solid line in the up- 
per panel indicate the observed AEz and the calculated AEz 
under the condition when Bn =0. The difference from the zero 
Bn line is the OHS at each Bz. As can be clearly observed, an 
abrupt decrease in AEz was observed at Bz =4.31 T owing 
to the bistable nature of NSP. In order to measure the degree 
of Bz compensation via Bn, we introduce the effective mag- 
netic field as experienced by e-spin; Beff(= Bz-Bn). By using 
the previously obtained values of g^ =+2.54 and gl =-0.37—, 
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we deduced the electronic splitting energy AEg = g^fiBBsSf, as 
shown in the middle panel of Fig. [He). In the region where 
fiz <4.31 T, the absolute value of AE^ nearly reduces to zero, 
and the Overhauser field fully compensates for the external 
field. With increasing B^, the magnitude of Z?n shows a clear 
reduction and lAiSgl increases abruptly. 

Here, we focus on the DCP of the PL (the lower 
panel of Fig. [TJc)). In this work, the DCP is defined as 
(r -/+)/(/"+/+) denotes the integrated PL intensity 
of the cr+(_) component). It is noteworthy that the DCP of 
PL shows a clear jump from ~0.6 to ~0.9; this transition 
synchronizes with the decrease in A^e- As mentioned above, 
the DCP of X^ is essentially determined solely by the e-spin 
polarization {^z), and it can be expressed as DCP=2(5'z). Ac- 
cordingly, a high (low) value of DCP indicates a small (large) 
reduction in e-spin polarization (i.e., e-spin relaxation). Since 
the e-spin relaxation rate depends on lAiSel (see Eq. [T]in the 
next section), the change in the X^ DCP presents the possibil- 
ity of a direct measurements of the electron and nuclear spin 
polarizations in a coupled e-n system. 

It is noteworthy that the DCP observed in the other charge 
states show different behaviors. The OHS observed in the 
X° and XX" PLs show changes similar to that observed in 
X^; this is one of the evidences that these PL lines originate 
from the same QDi^. In contrast, the tendencies of DCP are 
quite different for these other exciton complex peaks. The X" 
DCP stays constant (~0.7), thereby signifying independence 
from AEf., as shown in the lower panel of Fig. [He) (denoted 
by the gray symbols). This can be attributed to the contri- 
bution of the unpolarized X" supplied from the XX". XX" 
decays to X" by emitting cr+ and cr_ photons with identical 
probabilities, and therefore, XX'^ DCP is basically zero (not 
shown here). The DCP of X° is approximately calculated as 
[(« + n/4) - n/4]/[(« + n/4) + n/4] ~ 0.67, if QDs are ex- 
cited under the power at which n electron-hole pairs (X") are 
generated in each QD (in this case, XX°/X''=l/2 according 
to Poisson statics). Further studies in this direction requires a 
close examination of the fine structures of the exciton levels 
because of the electron-hole exchange interaction. 



B. Experimental estimation of correlation time and electron 
spin relaxation time 

We estimated the key quantities (t^ and Tjo) in the e-n spin 
dynamics from the experimental data. The e-spin relaxation 
rate under the effective magnetic field in frequency unit Qf^i- 
AE^Ifi) can generally expressed as 
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This equation represents a Lorentzian shape with a full width 
at the half maximum (FWHM) of in energy and an am- 
plitude of 1/tso. Here, Tc and Tso correspond to the e-n cor- 
relation time and the e-spin relaxation time at fig =0, respec- 
tively. As mentioned above, the e-spin relaxation rate appears 
directly in the reduction in X^ DCP, and these key quantities 
can be obtained from the fitting with the inverse of Eq. [T] 
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Electron energy splitting (|ieV) 

FIG. 2: DCPs as a function of the electronic energy splitting at the 
different external field. Symbols and colors indicate the experimental 
data and corresponding B^. Each of the data values was obtained 
from the time-resolved measurements of OHS and DCP. The absence 
of the data points around AE^ — 15 ^eV is attributed to the abrupt 
changes in OHS and DCP. The solid curve represents the Lorentzian 
fitting with a width of ~15 yueV and an offset of ~0.8. 



Figure |2] shows the DCPs as a function of A^e at differ- 
ent values of external field strength (2 T < Bz < 5 T). The 
data were obtained from the time-resolved measurements, as 
shown in a later section of the paper (Fig. [3jb)). As clearly 
shown, a definite dip is observed at AE^. ^ 0. The Lorentzian 
fitting (FWHM~15 fieV) with an offset (-0.8) depicted by a 
solid curve in the figure was able to reproduce the entire ex- 
perimental data. It should be noted that the X^ DCPs at the 
different Z?z values depict a unique curve, and this fact justifies 
the assumption that Tc and Tso are independent of Bz- 

First, we estimate the e-n correlation time Tc. By using the 
relation of FWHM= 2h/Tc, we evaluated the correlation time 
to be Tc ~80 ps. The obtained coincides with the values 
that Maletinsky et al2- (35 ps) and Braun et al? (50 ps) used 
in the calculations to reproduce their observations in single 
In(Ga)As QDs. To the best of our knowledge, an experimental 
estimation of this correlation time has not been reported thus 
far although Tq of several tens of picoseconds has been used in 
the data fitting. In addition, the value is in the same magnitude 
as the X^ decoherence time (~43 ps) that was measured by 
Fourier spectroscopy of the same InAlAs QDii. 

Secondly, we focus on the e-spin relaxation time under a 
zero effective magnetic field. This characteristic time can be 
expressed as Tjo - t^I [5op/5z (0) - 1 ji^, where 5op denotes 
the initial e-spin polarization injected into the QD ground 
state, and S z (0) is one when Beif =0. Here, tr denotes the 
recombination time, and it is found to be ~0.75 ns by other 
independent time-resolved measurements. From the DCP at 
A^e^O, 5'z(0) is evaluated to be ~0.6. In contrast, precise 
evaluation of 5 op is difficult because the e-spin polarization 
created with the cr_ excitation is lost partially during the en- 
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ergy relaxation process from the wetting layer to the QD 
ground state. For simpUcity, 5 op is replaced by the offset 
value of the DCP curve (~0.8) because the e-spin relaxation 
process is strongly suppressed in the region of large Z?eff ac- 
cording to Eq. [T] By using these values, we evaluated the 
e-spin relaxation time under A^e^O to be Tso~3rR. This value 
is in good agreement with the exciton spin relaxation time ob- 
tained in other measurements for resonant and nonresonant 
excitation&is— . 

Thus far, we were able to estimate Tc and t^q directly from 
the experimental data. Next, it is required to determine the 
factors influencing these characteristic times. Here, we con- 
sider the effects of the randomly fluctuating Overhauser field 
induced by the n-spin ensemble^. The random fluctuation 
of Bn (denoted by ABn) induces additional e-spin precession, 
and the memory of the e-spin polarization is lost from the ini- 
tial value. According to the traditional spin precession model, 
the longitudinal component of ABn (ABn,||) induces a loss in 
the transverse component of e-spin polarization {i.e., decoher- 
ence). On the other hand, the transverse component of ABn 
(ABn.j.) gives rise to a loss in the longitudinal e-spin polar- 
ization {i.e., relaxation). Assuming a random variable in NSP 
with a Gaussian distribution of a width y/N {N: the number 
of the nuclei in a QD), the fluctuation of Bn can be estimated 
as ABn = A/ y/Ng'^jdB {g'^' an isotropic electron g-factor. A: 
the hyperfine coupling constant). By using the typical values 
for a InAlAs QD {g". A, N) = (-0.37, 50 /leV, BxlO"*), we 
can roughly estimate the fluctuation of Bn to be ABn ~ 15 mT. 
The corresponding decoherence time (T^j ^) and the relaxation 
time Teno are given as the functions of ABn. 

To begin with, we compare ^ and the observed correla- 
tion time TiT^. The e-spin decoherence time induced by ABnji 

— and it is 

3 ns. The fact 
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can be expressed as 

estimated to have an approximate value of 
that the estimated Tj^^ is fairly longer than the observed Tc 
may indicate the presence of another scattering process which 
is responsible for the shorter decoherence time Tj. In our as- 
sumption, 1 /tc is given as a linear coupling of 1 / Jj „ and 1 /T^ 
and is dominantly determined by l/Tj. One of the plausible 
causes for l/T^ is the charge fluctuation in the QD region. 
Lai et al. have mentioned the effect of the e-spin tunneling 
between the QD and the n-doped layer in their charge-tunable 
QD structure-"*. Although our QD sample has no diode struc- 
ture and there is no interaction with the electrode, unlike their 
charge-tunable QD sample, similar phenomena such as car- 
rier tunneling between QD and its surroundings occur even in 
our sample; consequently various types of the charge states 
appear during the exposure time of the CCD detector (~0.1-1 
s), as shown in Fig.[TJa). The changes of the charge state may 
interrupt the e-spin precession, and it may affect the correla- 
tion time of an e-n spin system. Further Oresearch is required 
regarding the origin of l/T*. 

Next, we compare TenO and the observed e-spin relaxation 
time under B^s =0. Assuming an uniform electron wave 
function in the QD region, the e-spin relaxation time in- 
duced by ABn.j. under Beff=0, which is given as Teno - 

I^A^Tc (A/A^/!)^j , to be ~50 ns-^ Since the estimated Tgno 
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FIG. 3: (a) Current dynamics model of e-spin system including X* 
and the single hole states, and the corresponding PL polarizations, 
(b) Transient evolutions of DCP and OHS at =3.5 T. The exci- 
tation polarizations are depicted in the upper side of the panel. The 
solid curves are the calculated results in the coupled e-n system, and 
they were able to reproduce all the behaviors of the experimental 
results. Right inset indicates a schematic of the n-spin polarization 
(black solid curve) and depolarization (red dashed line) terms for an 
explanation of the transient OHS. 



is one order of magnitude longer than the observed t^q, we 
introduce the e-spin flip term (l/rgo) separately from the e- 
n flip-flop term (1/Teno), and the total e-spin relaxation rate 
obtained from the measurements is assumed to be I/Tjo = 
l/^enO + l/^eO- In our assumption, these two spin relaxation 
processes have the same Q.^ dependency as that represented by 
a common Lorentzian function. Note that the large reduction 
in the DCP shown in Fig. |2]cannot be reproduced without the 
introduction of I/t^o- 



In this section, the e-n correlation time and the e-spin relax- 
ation time at Beif =0, two key quantities in e-n spin dynamics, 
were evaluated from the experimental data. Since the esti- 
mated spin decoherence and relaxation time induced by the 
random fluctuation of Bn are of fairly longer duration than the 
values obtained from the measurements, other scattering pro- 
cesses are required to explain these durations. Although the 
source that decides and t^o has not yet been identified at 
the present stage, a direct evaluation provides the valuable in- 
formation for the modeling of e-n spin dynamics, as discussed 
in the next section. 
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C. Dynamics model of coupled electron-nuclear spin system 

Finally, we test the validity of the dynamics model of the 
coupled e-n spin system. The temporal evolution of the mean 
NSP ih) is described by the following rate equations^: 

d{h) 1 1 
-^^^m{S.)-So)-{h)]-—{h) (2) 

CZf 7 NF 1 ND 

where denotes the thermal e-spin polarization, l/T^Np 
and I/Tno denote the n-spin polarization and depolarization 
rates^i, respectively, and Q - /(/ + 1) / [5 (5 + 1)] denotes 
the momentum conversion coefficient from the e-spin to n- 
spin system. It is noteworthy that the n-spin polarization rate 
is also responsible for the e-spin relaxation, and it appears in 
the e-spin dynamics via l/rgno = A^/^nfo (1/7"nfo: the n-spin 
polarization rate at BefF =0). Although Eq. |2] with constant 
values of the averaged e-spin polarization (^z), has explained 
the observed OHS qualitatively in previous studies, the ac- 
tual (S'z) in the dynamics is expected to change along with the 
evolution of {I^. In our model calculation, {S^), which drags 
the randomly-oriented n-spin ensemble to the highly polar- 
ized state, is determined by the dynamics in the following four 
states, as shown in Fig.[3ja): with the spin-up/down elec- 
tron and «j: the populations of the corresponding states), 
and the spin-up/down single hole states (similarly denoted by 
«lt and and (^z) is given as (n-f - n[) / [2 (n-f + nj]. These 
four states are connected with the rates of the optical pumping 
with cr_ light, the radiative recombination (1/tr), and the spin 
flip of electron and hole (1/ts and l/rj,). 

The direct observation of the temporal evolutions of OHS 
and DCP can provide a better understanding of the e-n spin 
dynamics. Typical transients obtained from the target PL 
at Bz =3.5 T are shown in Fig.^b). In order to set the initial 
NSP to zero, the excitation polarization before the time of ori- 
gin was modulated between a+ and cr_, with a frequency of 
10 Hz. The temporal evolution of OHS {-A (I^)) is explained 
schematically by the difference between the n-spin polariza- 
tion (a black solid curve) and the depolarization (a red dashed 



line) rates, as shown in the right inset of Fig. Ob). After 
switching to cr_ excitation, OHS increases gradually in the re- 
gion of small difference and increases explosively around the 
peak of the n-spin polarization rate. Under this experimen- 
tal condition, the OHS jumps clearly to the saturated value 
within 3 s, and Bz compensation via Bn is achieved within the 
homogeneous broadening of the n-spin polarization rate. At 
the exact moment of the abrupt increase in the OHS, the DCP 
of the X^ PL drops suddenly from 0.8 to 0.7. The solid curves 
are the calculated results obtained from the abovementioned 
dynamics model. In the calculation, the key parameters (t^ 
and Tso) were in the same order as the experimentally evalu- 
ated values reported in a previous section of the paper The 
fact that the calculation could reproduce the observed DCP as 
well as OHS shows the validity of our dynamics model. 
IV. CONCLUSION 



In conclusion, we investigated the spin dynamics of the 
coupled electron-nuclear spin system in a single InAlAs QD. 
The DCP of X^ PL, which is basically determined by e- 
spin polarization, showed synchronized changes with the elec- 
tronic energy splitting, and this fact offers the possibiUty of 
NSP probing via X^ DCP in a QD structure. By taking advan- 
tage of this feature, the correlation time of the e-n spin system 
and the e-spin relaxation time, which play a crucial role in 
spin dynamics, were evaluated as ~80 ps and t^o ~3tr, 
respectively. The experimentally obtained Tc agrees well with 
the results obtained by calculations or the decoherence time of 
X^ PL; further, Tjo agrees well with the exciton spin relaxation 
time obtained from other experiments. Although the definite 
source of these key quantities were not identified at this stage, 
a direct estimation from measurements is very important in 
the characterization of the e-n spin dynamics. The spin dy- 
namics model used in this study successfully reproduces the 
observations of DCP as well as the OHS, and we believe that 
the model can significantly contribute to the understanding of 
e-n spin dynamics. 
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